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Abstract The temperature of fragmenting source in central heavy-ion collisions at Fermi energy is investigated by 


the isospin-dependent quantum molecular dynamics model in combination with the statistical decay model GEMINI. 


Five different nuclear thermometers are used to extract nuclear temperature. We find that the He and Li isotope 


temperature reaches a plateau at about 70-100 MeV/nucleon of beam energy. The slope temperature and the 


quadrupole fluctuation temperature give high values. The quantum slope temperature and the quantum quadrupole 


fluctuation temperature are more close to the He and Li isotope temperatures. 
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1 Introduction 


The nuclear matter may undergo a phase transition 
from the liquid ground-state to a gas of nucleons when 
the nuclei was heated. The nature of the interaction 
between nucleons is similar to the Van der Waals 
forces between molecules. The uncertainties of the 
nuclear equation of state (EOS) make the study of 
nuclear liquid-gas phase transition be important and 


meaningful”), 


Experimental signals of phase 
transitions have been observed by many experimental 
groups’*!, Many theoretical efforts to find signatures 
of phase transitions in heavy-ion collisions (HICs) 
have been performed a few years ago'”*!. HICs at 
intermediate energies offer an opportunity to heat 
nuclei, and the nuclear multifragmentation is long 
assimilated to the nuclear liquid-gas phase transition. 
Temperature is one of the most important 
degrees of freedom in describing phase transition. 
However, the nuclear matter created in HICs is a 
non-equilibrium, finite, and open system. It is more 


difficult to determine the temperature in nuclear matter 


than that in ordinary matter. Different from ordinary 


thermometer, nuclear thermometers are more 
complicated. Several nuclear thermometers have been 
proposed in the past days. These approaches can be 
divided into three families'”’” : (1) Kinetic approach. 
Kinetic approach is based on the assumption that 
energy spectra of particles obey the Maxwell- 
can be 


Boltzmann distribution. The temperature 


derived from the slopes of the kinetic energy 


[11,12] or the momentum fluctuations''?!, It has 


spectra 
been suggested that the slopes of light product spectra 
in nuclear reactions lead to very high “temperatures”. 
The temperatures extracted from Maxwell-Boltzmann 
kinetic approach probably reflect not only the thermal 
properties of the system, but also the collective 
energies coming from the dynamics of the nuclear 
collision. Recently, the Fermi-Dirac kinetic approach 
is considered to improve this type of thermometer” 4", 
(2) Population approaches. The underlying idea for 
this method is that the relative populations of the 
produced clusters or its excited states are assumed to 


obey a Boltzmann distribution. Population of excited 
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[16] [17] 


states ™ and double ratios of isotopic yields- ^- are two 
of the most often used methods. Employing the 
isotopic thermometer, Pochodzalla et al. had given the 
caloric curve which is taken as the evidence of the 
occurrence of a liquid- gas type phase transition"). (3) 
Thermal-energy approaches. The temperatures at 
freeze-out are obtained from the excitation energy 
which is extracted by measuring evaporation cascade 
varying 
neutron-to-proton ratio N/Z. An example is the isospin 
[18] 


from a thermalized source by 
thermometer 

In this paper, we attempt to compare five 
different nuclear thermometers. The central collisions 
of "°Xe+ "Sn at Fermi energy are simulated by the 
isospin-dependent quantum molecular dynamics 
(IQMD) model together with the statistical decay code 
GEMINI. The double ratio temperature, the slope 
temperature, the quantum slope temperature, the 
quadrupole fluctuation temperature, and the quantum 


quadrupole temperature are compared. 


2 IQMD transport model 


The IQMD model'’?”"! is based on the same principles 
as the quantum molecular dynamics model. With 
consideration of the different mean field potentials of 
proton and neutron, and the production of pion and 
kaon, the IQMD model has been widely and 
successfully used for the analysis of collective flows, 
stopping, pion and kaon multiplicities in HICs at 
energy below 2 GeV/nucleon”! The 
statistical code GEMINI was proposed by R. J. Charity 
in 1980s'*"!, Tt can be utilized to treat the decay of a 


incident 


compound nucleus in fusion reaction and excited 
fragments in HICs. GEMINI code was always applied 
to the transport model to statistically deexcite the hot 
fragments”, Recently, within the framework of the 
IQMD+GEMINI model, the odd-even effect in the 
yields of the final fragments has been well 
reproduced**!, Using the same model we have 
investigated the multiplicities and the kinetic energy 
spectra in central HICs, the simulations are in very 
good agreement with the experimental data. The slope 
temperature and the isotope temperature were also 
studied''*”4!, The transport with modified (the medium 
factor of 1- 0.2p / p,) nucleon-nucleon elastic cross 


section’! and soft EOS”®! are used to study the 


nuclear temperature. To compensate the fermionic 
feature the method of the phase-space density 
constraint in the constrained molecular dynamics 


model!” 


is applied. 
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Fig.1 Charge distributions of fragments produced in central 
collisions of '°Xe+'7°Sn at 50 MeV/nucleon. 


3 Results and discussion 


3.1 He and Li isotope temperature 


Figure 1 displays charge distribution of the central 
collision (b<1.6 fm) of 1X e4 Sn at 
50 MeV/nucleon. Both the simulation by IQMD with 
and without GEMINI code are displayed in the figure. 
The experimental data are shown by the symbol. 
Firstly, it is clearly seen that IQMD+GEMINI can 
reproduce the experimental data quite well. Secondly, 
with the help of GEMINI code, the yield of heavier 
fragments will decrease and that of light particles will 
increase, due to deexcitation of the hot fragments. The 
good agreement with experimental data, makes the 
study of the double ratio temperature be possible and 
believable. We 
temperature in which He and Li isotopes are used, 


introduce one of double ratio 


named Ty- It is defined by 


Y/Y 
Taa; =13.3MeV/ In| 2.2 5—4 |, (1) 
3 He “He 
where Y is the yield of the isotope. The T,,,,; can be 
calculated using this formula. 


3.2 The slope and the quantum slope temperature 


Figure 2 illustrates the calculated (open circle) and 
experimental (open star) kinetic energy spectra of free 
1Xe+!Sn at 


E\a=50 MeV/nucleon. The solid lines in each panel 


protons in central collisions of 


050513-2 


GUO Chenchen et al. / Nuclear Science and Techniques 24 (2013) 050513 


denote Maxwell-Boltzmann fit. The fit formula is 


Sel E) (2) 


slope 


expressed as 


Y(E) œ ( 


slope 


where Ep reflects the repulsive Coulomb forces”. 


Taking into account the Fermi nature of nucleon, the 


slope temperature was rewritten by Bauer”! 


A-A, 2 5m Tjee 
Tipe E -h pe ’ (3) 


where Ep is the Fermi energy. Tie names the 


slope 


quantum slope temperature. 
12Xe+!20Sn; b <1.6 fm 
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Fig.2 Center of mass kinetic energy spectra of free protons 
produced in central collisions of 29X04 Sn at Elap=50 
MeV/nucleon. 


3.3 Quadrupole fluctuation and the quantum 


quadrupole fluctuation temperature 


One of new nuclear temperature which is based on 
momentum fluctuations of detected particles was 
proposed in Ref.[17]. The momenta of particles were 
assumed to obey the Maxwell-Boltzmann distribution 


a 
1 7 e- = = P: } (4) 
2amtw, m 


fluct fluct 


f(p)= 


where m is the mass of the particle, px, py and p; are 
three momentum components of the particle. The 
variance o, can be obtained from the quadrupole 
Qy = Px- p, distribution through 


o: =|d°p(Q,) f (p) (5) 


Figure 3 shows the distribution of Q, for free 
protons. The form is similar to the experimental 


[13] 


results ~. Then, we can get the Tha from these 


uct 


formulae, the results are shown in Fig.4. More details 


will be discussed in the following section. The 
fluctuation temperature was improved as follows” ®: 
4E? 2r’, 
Tiuct = a + 15 Tas (6) 
where Er is the Fermi energy of nuclear matter and 
T! 


fluct 


temperature. 


is named as the quantum quadrupole 
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Fig.3 Distributions of quadrupole Q, = p- P; for free 


protons produced in central collisions of iaxe Sn at 
Eia» = 50 MeV/nucleon. 


3.4 Comparison between nuclear thermometers 


Within the IQMD +GEMINI model, we calculated the 
from "°Xe+"Sn 


collision at incident energies from 30 MeV/nucleon to 


nuclear temperatures central 


100 MeV/nucleon. Five nuclear thermometers are 
investigated in this work: the double ratio temperature 
Tuei, the slope temperature Tsiope, the quantum slope 


temperature T; fluctuation 


slope ? 


temperature Truc, and the quantum quadrupole 


the quadrupole 


temperature T; 


fluct 


. The results are shown in Fig.4. 
different 
thermometers give very different values of nuclear 


Firstly, it can be seen that 
temperature. The discrepancy between various nuclear 
temperature steadily grows with increasing incident 
energy. All temperatures exhibit an increase of value 
with increasing incident energy. Tą; reaches a 
plateau at about 70-100 MeV/nucleon of beam energy. 


and T, 


fluct 


Secondly, T 


slope 


are very close to each 
other, because both of these two approaches assumed 
energy distribution (Eq.(2)) and the momentum 
(Eq.(4)) are the 
Maxwell-Boltzmann type. These two temperatures are 


distribution classical 


the highest of all. It is because besides the thermal 
motion, they also contain the Fermi motion. Thirdly, 
Tie and Ty, 


slope fluct 


are both smaller than T. and 


slope 
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T, 


fluct * 


removed. But due to the different methods of 


This is because the Fermi motion of nucleon is 


elimination of the Fermi motion, the results are 
different. 


Temperature / MeV 


2Xe+2°Sn 


30 50 70 90 110 
E» / MeV-nucleon™* 


Fig.4 Nuclear temperatures extracted from various methods as 
a function of incident energies for central collisions of 
129¥ 64 1205, 


4 Conclusion 


Within the IQMD +GEMINI model, various nuclear 
thermometers have been investigated for the central 
collisions of '*°Xe+'?°Sn at incident energy range from 
30 to 100 MeV/nucleon. These nuclear thermometers 
are not exactly coincident with each other. The He and 
Li isotope temperature reaches a plateau at about 
70-100 MeV/nucleon of beam energy. The slope 
temperature and the quadrupole fluctuation 
temperature give higher values than others. This is 
because the Fermi motions of nucleons are not 
eliminated. The quantum slope temperature and the 
quantum quadrupole fluctuation temperature are more 
close to the He and Li isotope temperature because of 


the removed Fermi motion. 
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